Charge transport in a Tomonaga-Luttinger liquid: effects of pumping and bias 



Amit Agarwal and Diptiman Sen 
Center for High Energy Physics, Indian Institute of Science, Bangalore 560012, India 

(Dated: February 6, 2008) 

We study the current produced in a Tomonaga-Luttinger liquid by an applied bias and by weak, 
point-like impurity potentials which are oscillating in time. We use bosonization to perturbatively 
calculate the current up to second order in the impurity potentials. In the regime of small bias and 
low pumping frequency, both the DC and AC components of the current have power law dependences 
on the bias and pumping frequencies with an exponent 2K — 1 for spinless electrons, where K is 
the interaction parameter. For K < 1/2, the current grows large for special values of the bias. For 
non- interacting electrons with K = 1, our results agree with those obtained using Floquet scattering 
theory for Dirac fermions. We also discuss the cases of extended impurities and of spin- 1/2 electrons. 
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I. INTRODUCTION 

The conductance of electrons in a quantum wire has 
been studied extensively in recent years both theoreti- 
cally [H, i] and experimentally [3, i, i, i, 0, i- For a 
wire in which only one channel is available to the elec- 
trons and the transport is ballistic (i.e., there are no im- 
purities inside the wire, and there is no scattering from 
phonons or from the contacts between the wire and its 
leads), the conductance is given by G = 2e^//i for in- 
finitesimal bias. However, if there is an impurity inside 
the wire which scatters the electrons, then the conduc- 
tance is reduced. For a (5-function impurity with strength 
[/, we obtain G — (2e^//i) (1 — U'^/vp), to lowest order 
in U, where is the Fermi velocity of the electrons. 
In the presence of interactions between the electrons, the 
impurity strength U effectively becomes a function of the 
length scale through a renormalization group (RG) equa- 
tion [il, • The RG flow has to be cut-off at the smallest 
of the three length scales of the system, namely, the wire 
length, the thermal length which is inversely proportional 
to the temperature, and a length which is inversely pro- 
portional to the bias voltage Vbias- If the latter length 
scale is the smallest of the three, then the effective value 
of the impurity strength is given by J/VJ^^^^i where K 
is a parameter related to the strength of the interac- 
tions between the electrons as we will see later. Hence 
the correction to the conductance due to the combined 
effect of the impurity and the interactions is given by 
AG - C/^K^^"^ 

^ * bias 

The phenomenon of charge pumping and rectification 
by oscillating potentials applied to certain points in a 
system has also been studied theoretically [11-38] and 
experimentally [H, HO, El, 113, El . For the case of non- 
interacting electrons, theoretical studies have used adi- 
abatic scattering theory 0, [l^ , Floquet scattering 
theory |23|, i2^] , variations of the non-equilibrium Green 
function formalism [2^, [2^, , and the equation of mo- 
tion approach fssj. The case of interacting electrons has 
also been studied, using a RG method for weak interac- 
tions [31 , and the method of bosonization for arbitrary 
interactions [45-55]. The analytical methods used in the 



two cases typically treat the electrons in quite different 
ways, with a non-relativistic Schrodinger equation or a 
tight-binding model on a lattice being used in the non- 
interacting case, and a massless Dirac model followed by 
bosonization being used in the interacting case. 

A clear comparison between the non-interacting and 
interacting cases does not seem to have been made be- 
fore. We plan to fill this gap in this paper and will study 
the effects of electron-electron interactions (of arbitrary 
repulsive strength) on the DC and AC components of 
the current in a system with a bias and time-dependent 
impurities. 

In Sec. II, we discuss a massless Dirac model for non- 
interacting electrons in the presence of several point-like 
impurities. We use Floquet scattering theory [2I, [24| to 
study the pumped and bias current in this model. In 
Sec. Ill, we first review the calculation of the backscat- 
tered current. Using bosonization [H, 113, HI], we then 
compute the DC and AC components of the current up to 
second order in the impurity potentials [H,!!^]- We show 
that these reduce to the results obtained using Floquet 
scattering theory for the non-interacting case. The cases 
of extended impurities [13] and of spin- 1/2 electrons are 
discussed briefly. We summarize our results in Sec. IV. 
We will consider infinitely long wires and zero temper- 
ature throughout the paper; hence the relevant energy 
scales in the problem are set only by the bias and the 
pumping frequency. 



II. MASSLESS DIRAC MODEL WITH 
IMPURITIES 



In this section we consider a spinless and massless 
Dirac fermion with no interactions between the fermions, 
and use the Floquet scattering theory to compute the to- 
tal current. The Hamiltonian in the presence of several 
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point-like impurities is given by H = Hq + Himp, where 



Ho 



Himp ^ dx 22 ^{x-Xp) Up{t) %l)\x)il;{x), 

J p 

i;{x) = ^h{x) e^'"^ + ^bLix) e-^'=-^ (1) 



where ipL and ipR are the fermionic field operators of the 
left and right moving electrons, is the Fermi velocity, 
and kp is the Fermi wavenumber which originates from 
some underlying microscopic model. For instance, one 
may have a system of non-relativistic electrons with a 
Fermi energy Ep = kp/{2'm) and vp = kp/m. (We are 
setting Planck's constant h equal to unity) . The Hamilto- 
nian Hq is obtained by linearizing the dispersion around 
the two Fermi points given by k — ±kp. Himp arises 
from the time-dependent impurities which have strengths 
Up(t)] this Hamiltonian couples left and right moving 
fields since 

(2) 

We will assume that 



Up{t) = Up cos{ujt + 4>p) , 



(3) 



i.e., all impurities vary harmonically in time with the 
same frequency lu. We will now use Floquet scattering 
theory and carry out a perturbative expansion in the di- 
mensionless quantities Up/vp. 

The equations of motion in the presence of a single S- 
bmction impurity 5{x — Xp)Up cos(ajt (j)p) is as follows: 



IVp 



dt ' " dx 
-- 5{x-xp) UpCOs[ujt + (j)p) ii^R + ijLe'''^'""'''), 



IVp 



dt dx 
-- 5{x-xp) UpCOs{ujt + (j)p) (V'L + V'fle'^''''"''')- 



(4) 

If we define the linear combinations 1^+ = T/'^e*'^^^'' + 



~dt 



dip+ „ 
ivp — — = 0, 
ox 



dip dip 
^ + i'^F ^ — — 26{x — Xp)UpCOs{ujt + (j)p) tp^. 

(5) 

By integrating over a little region from Xp — e to Xp + e, 
we find that ipj^ is continuous at the point x = Xp^ while 
-0- has a discontinuity given by 

ivp [■(/;_ (xp-fe) — V-'-(a;p—e)] = 2UpC0s{Ljt-\-(j)p) i)^{xp). 

(6) 



We would like to note here that it is necessary to retain 
the terms V'ljV'i?, + V'lV'l in Eq. Q in order to have con- 
tinuity of Tp+ . In some papers, these terms are not taken 
into consideration. One then runs into the mathemati- 
cal peculiarity that ip+ and V'- are both discontinuous 
at X = Xp, and the discontinuity is taken to be propor- 
tional to their values at that point; but those values are 
actually ill-defined due to the discontinuity. 

We can now solve Eqs. (jH) along with the bound- 
ary conditions in Eq. For a single (5- function im- 
purity oscillating with frequency a; at x = Xp, let us 
consider a wave coming from the left {x < Xp) with 
energy Eq and unit amplitude. Note that we are mea- 
suring energies with respect to a Fermi energy, so that 
-©0 = corresponds to a fermion at the Fermi energy. 
Due to the oscillating impurity potential, the wave will 
be reflected back to the left with energy En = i?o + nu) 
and amplitude SLL{En, Eq), or transmitted to the right 
(a; > Xp) with energy En and amplitude SuLiEn, Eq), 
where n — 0, ±1, ±2, • • • defines the Floquet side bands 
[2^ . Note that since we are considering a Dirac fermion, 
there is no upper or lower bound to the energy En, and 
the velocity uj;^ is independent of the energy. (This is 
unlike the case of a non-relativistic fermion or a fermion 
on a lattice where there is a lower or upper bound to the 
energy, and the velocity is a function of the energy) . To 
be explicit, the wave function is given by 

= gi(kox~Eot} f^^, x<Xp , 

= ^ SRL{En,EQ) e'ik^x-E^t) a;>xp, 



i^L ^Y. 5ll(S«,So) e'(-'="^-^"*) for x < Xp 



= for X > 



jp , 



(7) 



where fc„ = En/vp. Similarly, we can consider a wave 
coming from the right with energy Eq and unit ampli- 
tude; it will be reflected back to the right with amplitude 
SR.R{En, Eq) Or transmitted to the left with amplitude 
SLR{En, Eq). Let us simplify the notation by defining 

rL.n — SLL{En, Eq), tL,n — SlrXEu, Eq), 

tR.n = SRL{En, Eq), rR,n — SRR{En, Eq). (8) 

Due to unitarity, we have the relations 

n 

E + i^i.-n = 1- (9) 

n 

The different Floquet scattering amplitudes rQ,_„ and 
ta^n can be found by using the boundary conditions in 
Eq. We will consider the case of several impuri- 

ties labeled by the index p as in Eq. U). To simphfy 
our calculations, we will assume that uj{xp — Xr)/vp and 
Eq{xp — Xr)/vp are small for all pairs of impurities p and 
r; the first condition corresponds to the adiabatic limit. 
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while the second condition impHes that we are only con- 
sidering states close to the Fermi energy. Keeping terms 
only up to first order in Up/vp^ we find that only the first 
Floquet side bands are excited, and 



tL. 



I 

i 

2v^ 



p 



~i4>p 



TLA 



rR,i 



rR,-i 



2vF 



^i{2kpXp — cl>p) 



P 



2vF 



E^. 



^i( — 2kpXp — cj>p) 



^i(~2kFXp+(t>p) 



(10) 



We also find that the unitarity relations in Eq. ^ are 
satisfied up to second order in Up/vp, and therefore tL,o 
and tj^ Q are given by 



\tL,o\' = l-\rRM^-\rB..-i\'-\tLA\'-\tL^^i\^ 



\tR,of = 1 



\tR,lf - \tR,-l\ 



(11) 



to that order in Up/vp- Note that the amplitudes given 
in Eqs. (jlOlllip are all independent of Eq under the ap- 
proximations that we have made. 

The depart of the current in, say, the right lead is 
given by [23l | 



iR,dc 



/OO 7 pi 
E i\^H,n\'{fR{Eo)-fR{E^)} 



+ \tR,n\^ {fUEo) - fR{E.n)} 



(12) 



where q is the charge of the electron, and fa{E) = 
l/[e(^~''°^/'''^-^ -I- 1] is the Fermi function in the lead a. 
At zero temperature, fa{E) = d{^a — E), where 9{E) = 1 
for E > and for E < 0. Note that in a non-relativistic 
system, an expression like (IT^ would contain ratios of ve- 
locities Vn/vQ. Since we are considering a massless Dirac 
fermion here, the velocity is independent of the energy, 
and Vn/vQ = 1 for all n. 

Let us define a frequency in terms of the bias voltage, 
luq = qVbias = ^J'R — IJ^L- We will assume the bias to be 
small so that ujQ{xr — Xp)/vp is small for all values of p 
and r. Using Eqs. (|10j|lip . we find, up to second order 



in Up/vp, that 
qujQ 

271 

qw 
2^ 



lR,dc = - ^ [\tR.of + \tR,lf + \tR,-lf] 

[\tR.i\^-\tR.-i\' + \rR,i\'-\rn,^i\ 



qwo qujQ 



2n Anvp 



E^P 



2 E] UpUr C0s{2kpXrp) COs{(j)rp)] 



quj 



E] UpUr sm{2kpXrp)srcL{(j)rp)^ (13) 



where 



and 



. Eq. (Uni) shows 



the effects of a bias (wq) and harmonically oscillating 
potentials (w). For the pure pumping case with wq = 0, 
Eq. (fT3|l agrees with the results presented in Ref. [ssl ]: 
note that the pumped current depends on sin((/)rp). 

It is interesting to note that the first term is just the 
ballistic conductance of a clean wire multiplied by the 
bias, the second term is a correction to the clean case be- 
cause of the presence of impurities, and the third term is 
the pumped current. In the non-interacting case, the bias 
component and the pumped component separate out, but 
for the interacting case, the current involves powers of 

Wo ± OJ. 



III. BOSONIZATION CALCULATION OF 
BACKSCATTERED CURRENT 

A. Backscattering current operator 

We now compute the current in a system of interact- 
ing electrons using the backscattering current operator 
introduced in Refs. [H [H [Ulil . 

Let us take the impurity potentials to be absent at 
time t = — oo; then they are gradually switched on. At 
the initial time, Hq commutes with the number operators 
of the left moving and right moving fermions, and 
Nfi respectively. In the absence of any impurity poten- 
tials, all the right movers originate in the left reservoir 
which is maintained at the chemical potential /ul, and 
all the left movers originate in the right reservoir main- 
tained at the chemical potential fiR. Hence, the system 
is initially described in the grand canonical ensemble by 
the chemical potentials /il and fiR which are the coeffi- 
cients of the number operators and Nji respectively. 
We will work in the interaction representation, taking the 
chemical potentials to be part of the interaction. This in- 
troduces time dependences into the fermionic operators 
tpL — ^ V'lc*'^^* and and ipR i/jfle*^"*. The operators 
tp^fpR and Tp^nipL appearing in in Himp (see Eqs. ^ and 
^) therefore pick up factors of e'''*'^°*. 

If there were no impurities, there would be a cur- 
rent flowing to the left given by Jq = g^Vbias /(27r) = 
qujo/ (2Tr). In the presence of impurities, some of this 
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current is backscattered to the right. The total current 
flowing to the right is given by / = —Iq + I^s, where I^s 
is the correction to the current due to backscattering by 
the impurities. The backscattered current is defined as 

lbs (t) q = -iq [Nr , Htmp ] 

= iqY. Up{t) [ Ai^R e'('^°*+''=-"'') - H.c. ]. 
p 

(14) 

The backscattered current at any time t is given by 

{hs) = (0| 5(-oo;t) hS) S{t--^) |0), (15) 

where |0) denotes the initial state at t ^ — oo, and S is 
the scattering matrix arising due to the impurities, 

S{t]—co) = S'^(— oo;t) = T cxp[— i / dt' Himp(t')]. 

J — oo 

(16) 

We define a backscattering operator 

In terms of this, Ibs{t) ~ ^Q'^p[Bp{t) — H.c], while 

5(t,-oo) = 1 - iJ^pI-^dt'lIBpit') + H.c] to first or- 
der in Up. Thus 

Its ^ qY. f [Bp{i')Bl{t)-Bl{t')Br{t)] 

+ H.C |0) (18) 

to second order in U^,. 



B. Bosonization 

In one dimension, it is known that a large class of 
fermion systems which are gapless and have a low-energy 
dispersion which is linear can be written in terms of gap- 
less bosonic systems [HBjIs^iIH^- These systems are called 
Tomonaga-Luttinger liquids; for spinless fermions, they 
are characterized by an interaction parameter K and a 
velocity v. Non-interacting fermions have K ^ 1 and 
V = vp, while K < \ corresponds to repulsive interac- 
tions between the fermions. To be specific, consider a 
system with short-range density-density interactions of 
the form 



Hi, 



dxdy p{x)V{x - y)p{y), (19) 



where V{x) is a real and even function of x, and the den- 
sity p = ^^tp is given in Eq. We can write Eq. ITOl) 
in a simple way if V{x) is so short ranged that the argu- 
ments X and y of the two density fields can be set equal 



to each other wherever possible. Using the anticommu- 
tation relations between the fermion fields, we obtain 



Hi. 



92 



dx IptlpRll^lll^L, 



(20) 



where 32 is related to the Fourier transform of V{x) 
as g2 = V{0) — V{2kp). Defining a parameter 7 = 
g2 / {'2,t:v f) , we have the relations 



K 



1-7 



1/2 



1 + 7. 



(21) 



by 



In the absence of impurities, the bosonic action is given 



S = 



dtdx 



1 

2^ \dt 



V f d(j) 
2 Ife 



(22) 



Bilinears in fermion operators can be written in terms of 
bosons and Klein operators [HI, [13, Ull , such as 



2na 



~t - -i2y/7Kd, 

VlVR e ^ 



(23) 



where rjR and rjL are the Klein operators, and a is a 
short distance cut-off. We then obtain the ground state 
expectation value of products of four fermion operators 
as in Eq. p^ . namely. 



IplliXp, t')lpL{Xp, t')tpl{Xr,t)lpR{Xr,t) |0) 



(27r)2 [{xp - Xrf - {v(t' ~t)~ ia) 



2]K 



(24) 



for all values of K. 

For the non- interacting case with K = 1, we can eval- 
uate the above ground state expectation value directly 
without using bosonization. We use the second quan- 
tized expressions for the fermion fields. 



-tpR 



— 00 
00 



dk 

2^ 
dk 
2^ 



a-Rk e 



ik{x — VFt) 



a-Lk e 



ik{~x—vp-t) 



(25) 



where the creation and annihilation operators satisfy the 

anticommutation relations {dRkjORk'} ~ {'^Lk, o,]^k'} ~ 
2TT6{k — fc'). The ground state |0) is annihilated by 
ciB.kjCiLk for k > and by a]^j,,a^j, for fc < 0. We then 
find that the ground state expectation value agrees with 
the result given in Eq. ([M)) for if = 1 and v — vf- 

In general, the backscattered current has two parts: 
one independent of time which we call J^c, and the other 
varying with time, with frequency 2iu to second order in 
Up, which we call lac- lac does not contribute to any 
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charge transfer as its average over a cycle is zero. In the 
next few subsections, we calculate the expectation value 
of the backscattered current for various cases and study 
them in different limits. To simplify our calculations, we 
again assume that ojXrp/v and ujQXrp/v are small and that 
> 0. It will be convenient to define the combinations 



(JJ+ — ojQ + uj, and lo- — ojq 



(26) 



a) as Up{L) = {L/aY~^Up. In our case, the length 
scale L is set by or v/\uj-\. The effective impu- 

rity strength Up{L) therefore increases as {v /\Lo±\Y~'^Up 
for Ljj^ or LD- — > 0, and the correction AG grows as 
[Up{L)]'^ ^ Up\uj±\'^^~^ . This divergence must be cut off 
when AG becomes of order 1, in units of / {2tt). Restor- 
ing the appropriate dimensionful quantities, we see that 
the above RG analysis and perturbative expansion are 
vahd as long as Up/v « {a\uj±\/vy~^ . 



C. Single impurity 

This case has been discussed in Ref. [48|; we repro- 
duce the results here for the sake of completeness. Some 
details of the calculations are provided in the Appendix. 



jpp ^ lUp 



(- 



2K-2 



X [sgn{uj+) \uj+ 



\2K-1 



+ sgn{uj-) \u!-\ 



2K-1] 



(27) 



(?) 



2K-2 



8Trv^T{2K)cos{TrK) w. 

X cos{2ujt + 2(f)p + sgn{uj+)TTK) 
+ sgn{LU_) 

X cos{2LLjt + 2(t)p — sgn{u!-)TTK)], 



(28) 



where sgn{n) = 1 if > 0, if = and -1 if < 0. 
In Eqs. (|27ll28p . we note that the currents become large 
in the limit loq ii-j if K < 1/2. Hence the perturba- 
tive expansion in powers of Up breaks down when loq is 
close to ztw [il]. The region of validity of the perturba- 
tive expansion can be estimated using a RG analysis as 
discussed below. 

Eqs. (|27ll28p imply that for the pure pumping case 
with Wo = 0, Ij^s.dc = ^bs.ac = 0- For a single impurity, 
therefore, charge pumping does not occur, whether or 
not there are interactions between the electrons. For the 
pure bias case with a; = and </)p = 0, we have /^f -I- 
^bsac ^p^V^~^- Thus the backscattering correction to 
the conductance given by -Ibs,dc/Vbias = -qIbs,dc/uJo is 
proportional to C/p V^^f,"^ 

In the presence of both bias and pumping, the correc- 
tion to the differential conductance AG — —qdIbs,dc/duJo 
grows large as Up\uj±\'^^~'^ for uj^ or aj_ — > 0. This is 
consistent with results based on RG calculations 0, • 
Namely, the presence of interactions between the elec- 
trons effectively makes the impurity strength Up a func- 
tion of the length scale; this is described by the RG equa- 
tion dUp/d\nL = (1 — K)Up, to first order in Up{L). 
Hence the value of Up{L) at a length scale L is related 
to its value Up defined at a microscopic length scale (say, 



D. Several impurities 

We now consider the case of several impurities located 
at Xp with the phases of the oscillating potentials being 
(j)p. We again define x^p and (prp as in Eq. (|f 3p . The 
backscattered current can be written as Ihs — ^ts "I" 



J2p<r^bs- The dc and ac parts of /^f are given in the 
previous subsection. Next, we find that 



qUpUr 



2K-2 



At:v'^T{2K) 
X [sgn{ujj^) \loj^\'^^^^ cos{2kpXrp + (t^rp) 
-\- sgn{ui-) \ui-\'^^~'^ C0s(2kFXrp — (t>r 



(29) 



qUpU, 



€,ac = 7^::57v^n&7rFT (^) cos(2fc;^a;,p) 



ATrv^T{2K) co&{t:K) 
X [S5n(c.+) 

X cos(2ujt + 4>p + (f)r + sgn{u)+)TTK) 
+ sgn{uj-) 

X cos{2u!t + <j)p + — sgn{LU_)TTK)]. 



(30) 



For the pure pumping case with loo — 0, we see that 
^bl.dc ^ ^^■^"^ sm{2kFXrp) sm{(prp), while I^^^ = 0. Eq. 
(I^^l differs from the results given in Ref. due to the 
terms involving 2kFXrp- 

We note that the currents given in Eqs. (|27ll28p and 
(|29ll30p all reverse sign if we change — *■ — (^o ^^nd Xp — > 
—Xp for all p. This is a natural consequence of parity 
reversal, i.e., interchange of left and right. 

The dc parts given in Eqs. ((27|) and ([29|) can be 
combined to give a total current /ftg = X]p ^bs dc + 



jpr 
p<r bs.dc'> 



^bs.dc 



8TTv^r{2K) 



2K-2 



X [sgnicj+) |c.+ |2^-i I J2Upe 



i{2kFXp+4>p) |2 



-I- sgn{uj^ 



\2K-1 



p 



JJ^gi(2kFXp~4>p) |2 



(31) 
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The above expression suggests that the current will be 
maximized if either 2kpXp + (j)p or 2kpXp — 4>p has the 
same value for all p. This means that the potentials 
in Eq. ([3]) should be of the form UpCos{u)t — 2kpXp) 
or UpCos{Lut + 2kFXp); this describes a potential wave 
traveling to the right or to the left. Such a wave has 
been studied extensively for the case of non-interacting 
electrons; see Refs. [M, HI HI iH, IH, f33, 34, 38] and 

An unusual phenomenon occurs if the interactions are 
sufficiently strong, i.e., if i^T < 1/2. If there is no bias, 
the DC part of the current generally goes as u)^^~^ which 
increases as uj decreases. However, it is clear that if lo 
was exactly zero (time- independent impurities), then the 
current would also be zero. These two statements imply 
that the current must be a non-monotonic function of oj, 
and must have at least one maximum at some value of 
UJ. Finding the location of the maximum requires us to 
go beyond the lowest order perturbative results of this 
paper. 




FIG. 1: (Color online) DC part of the backscattered cur- 
rent as a function of the bias uiq for several impurities, 
when uj± are small. The red (dot), magenta (dash), 
blue (dash-dot) and black (solid) lines show the results 
for K = 1/4,1/2,3/4 and 1 respectively. We have taken 

Figure 1 shows the dc part of the backscattered cur- 
rent as a function of the applied bias, for a fixed non- 
zero value of the pumping frequency assuming that 
LJ and wq a-re small. We have used the expression in 
Eq. ([31]) to plot the value of Ibs,dc{^Q)/Ibs,dc{^Q = 0) 
as a function of wq/w, for four different values of the 
parameter K = 1/4,1/2,3/4 and 1, taking the ratio 
I Ep C/pe'(2fef |2/| [|^gi(2fcp:cp-0p) |2 = 2 : 1 as an 

example. For K — l/A^ the current diverges at loq ~ ±a; 
as mentioned above. We also note the linear and piece- 
wise constant dependences of the current on for K = 1 



and 1/2 respectively; this is discussed in Subsec. III. E 
below. 

If we relax the assumptions that Loxrp/v and uoQXrp/v 
are small, then the exact expressions (up to second order 
in the impurity potentials) for the ac and dc component 
of the backscattered current are given by 



jpr 

bs.dc 



hrp I 



1/2-K 



X [S9n{uj+) \uJ+\^^^/'^JK-l/2{\^+Xrp\/v) 



X cos(2A;Fa;rp + 4>rp) 
-\- sgn{uj-) \uj- 

X C0s{2kpXrp 



"^-'Z' JK-l/2{\^~Xrp\/v) 
' 0rp)], 



(32) 



jpr 

bs,ac 



qUpUr^/n 



ATrv^T{K) cos(7rii:) 



TlK) \ vJ 



2K-2 



1/2-K 



COs{2kFXrp) 



X [sgn{uj+) JK_i/2{\^+Xrp\/v) 
X cos(2wt + (f>p + (j)r + sgn{uj+)'KK) 
-f sgn{uj^) \uj^\^^'^/^ JK-i/2{\^-Xrp\/v) 
X cos(2a;t + (f>p + (pr — sgn{u}-)'KK) 

+ J^^^_j^(\^^xrp\/v) 

- J,/,_K{\io.Xrp\/v)} 
X sin{2ujt + (f>p + (l)r)]. (33) 

The Bessel function J is discussed in the Appendix; us- 
ing a power series expansion given there, we can show 
that Eqs. (|32ll33p reduce to Eqs. (|29ll30p in the limit 
uj±Xrp/v 0. We note that the expressions in Eqs. ([771 - 
do not change if we relax the assumptions that oj and 
luq are small. 

Figure 2 shows the dc part of the backscattered current 
as a function of the applied bias for the case of two im- 
purities, labeled 1 and 2, taking Ui — U2, 2kpxi2 ~ 7r/2, 
(j)i2 — — 7r/4, and UJX12/V = 1; thus lo and ojq are not 
small, in contrast to the case shown in Fig. 1. We 
have used the expressions in Eqs. ([77]) and ([5^ to 
plot the value of Ibs,dci^o) / hs.dci^a = 0) as a func- 
tion of loq/uj, for four different values of the parameter 
K = 1/4,1/2,3/4 and 1. We see some oscillations in 
Fig. 2 due to the appearance of the Bessel functions in 
Eq. For K = 1/4, we again see divergences at 

UJQ = ±UJ. 

E. K = 1 and 1/2 

We now discuss the special cases K = 1 and 1/2 where 
the expressions for some parts of the currents simplify 
considerably. 
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FIG. 2: (Color online) DC part of the backscattered current 
as a function of the bias ujo for two impurities, when uj± are 
not small. The red (dot), magenta (dash), blue (dash-dot) 
and black (solid) lines show the results for K = 1/4, 1/2, 3/4 
and 1 respectively. We have taken Ui — U2, 2kFXi2 = vr/2, 
012 = — 7r/4, and LOX12/V — 1. 



and (I29II30I1 . We find that 

-I- sgn{uj-)], 
[sgn{uj^) cos{2kpXrp + 4>rp) 



jpp 

bs,dc 



rPr 

bs,dc 



I; 



jpr 



[sgn{uj+) 

Hwav 

qUpUr 



Anav 



pp _ 



^ ^ -[{sgn{uj 



Swav 



Anav 



+ sgn{ujJ) cos{2kFXrp — (prp)], 

- sgn{u!^)) cos{2ujt + 2(/)p) 



- ln| — I sin{2ujt + 2(j)p 



[(sgn(w+) -I- sgn{uj^)) cos(2ajt + (fip + (pr) 



- ln| — I sm{2ujt + <f)p + (})r)]. 



(38) 



Thus the DC parts of the currents do not depend on the 
precise values of w and Wq if they are unequal, and they 
have a finite discontinuity when uj crosses icjo. 

To conclude, we see that the dc parts of the currents 
are linear functions oiuJojUj for K = \, and are piecewise 
constant functions of 0^0,0; for K = 1/2. 



For non-interacting fermions with = 1, we find from 
Eqs. I|27ll28p that in the single impurity case, 



tPP _ ^ P , , 



qU^ 



Iliac ^ i^^O C0s(2^t-f 20p). 



(34) 



The total current is given by / = -/q + I^sdc + ^bsac^ 



I = 



qtJo 
2tt 



- 1 



Up cos(ijji -I- 4>p) 



VF 



(35) 



This is consistent with the fact that the transmission 
probability across a static point-like barrier of height U 
is 1 — {U /vfY' up to order . For the case of several 
impurities, we find from Eqs. (I29j|30p that 



jpr 
bs.dc 



jpr _ 



— — 5- uja cos[2kFXrp 



cos 



'rp)j 

UJ sm{2kFXrp) sin((/)r 



Lrp) aiLiyi^rp 

^^^^2 "^0 cos{2kFXrp) cos(2cjt 



(36) 
(37) 



Note that the dc part of the current is given by a linear 
combination of the pure bias part and the pure pumping 
part, and it agrees with the expression given in Eq. (|13p . 

For K = 1/2, we can obtain the different parts of the 
currents by taking the limit K 1/2 in Eqs. (|27ll28p 



F. Extended impurities 

The analysis in Subsec. III. D can be readily gener- 
alized to the case where there is an extended region of 
oscillating potentials [H^l • Let us replace the discrete set 
of potentials given in Eq. ([3]) by an oscillating potential 
of the following form 



C/(t) 



dx U{x) cos[ujt + (t>{x)] 



(39) 



We then see from Eq. (|3ip that the dc part of the 
backscattered current is given by 



IbsAc — 



(?) 



2K-2 



Trv^T{2K) 

X [sgn{Lu+) \lo+\^^-^\ J dx U{x) e^[2fcF^+'/>W] p 

+ sgn(w_) \uJ-\'^^-^\ J dx U{x) e»[2fc^^-'^(^)l p]. 

(40) 

to second order in U (x) . For the pure pumping case with 
luq = 0, wc find that 



bs.dc 



(- 

V V 



2K-2 



.2K-1 



At:v'^T(2K) 

^ j j dxdx' U{x)U{x') sin[2fcF(x- x')] 

X sm[<f>{x) ~ (p{x')]. (41) 

Eq. (liTI) implies that the charge pumped per cycle, 
AQ — {2n/cj)Ibs,dc, scales as w^^"^; for K < 1, this 
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grows large in the adiabatic limit w ^ 0. In this limit, we 
saw earlier that the effective length-dependent impurity 
strength diverges at small energy scales, which implies 
that the impurity presents a very large barrier to the 
electrons and the transmission coefficient is very small. 
In this limit, it has been argued in Refs. [3,143] that the 
pumped charge AQ is quantized to be an integer multiple 
of q. 



G. Spin-1/2 electrons 

For spin-1/2 electrons in one dimension, the phe- 
nomenon of spin-charge separation occurs if there are 
interactions between the electrons. The spin and charge 
degrees of freedom can be separately bosonized [s^ [s^] . 
The two bosonic theories are characterized by the pa- 
rameters {Ks^Vs) and {Kc,Vc) respectively. For a system 
with SU{2) rotational invariance. Kg = 1. The ground 
state expectation value in Eq. ([M]) then takes the form 

1 



[{xp-xrf ~ {vs{t' ^t)-iay]^/^ 

1 

" [{Xp~Xry ~~ Mt' - t) - ia)^]K^/^ 



, (42) 



where a is the spin label. The appearance of two 

different velocities, Vs and Vc, and two different expo- 
nents, 1/2 and Kc/2, in Eq. makes the expressions 
for the backscattered current rather complicated. How- 
ever, we can find the power law of the dependence of 
the currents on the frequencies by a simple scaling argu- 
ment. With the approximations made earlier, ujXrp/vs^c 
and ujQXj-p/vs i. — > 0, we see that the time dependence has 



changed from l/{t' — tf in Eq 



to l/{t'-t) 



Eq. ([^^ . The dependences of the backscattered currents 
on the frequencies therefore change from |a;o±a;p^~^ in 
the spinless case to |a;o±w|^'= in the spin-1/2 case. Since 
Kc is positive in general, the current no longer diverges 
as ujQ ±07. 



IV. DISCUSSION 

We have considered the effects of a bias and a number 
of weak and harmonically oscillating potentials on charge 
transport in a Tomonaga-Luttinger liquid. We have com- 
puted the backscattered current to second order in the 
amplitudes of the potentials. For most of our results, we 
have assumed the oscillation frequency and the bias to 
be small, but we have relaxed that assumption in Eqs. 
(|32j|33p . For our assumption of a Dirac fermion with a 
linear dispersion to be valid for an experimentally real- 
izable system, we must of course assume that uj and wq 
are small compared to the band width of the electrons. 

We find that the backscattered current is maximized 
for a traveling potential wave in which the positions and 



phases of the oscillating potentials are related in a linear 
way. For spinless electrons, if the interactions are suffi- 
ciently repulsive with K < 1/2, the backscattered cur- 
rent diverges for special values of the bias, namely, for 
ujQ — > icj. For any repulsive interaction, with K < 1, 
the correction to the differential conductance diverges 
for Wo — > iij-". Finally, we have pointed out a pecu- 
liarity which arises when several impurities are present 
and K < 1/2; namely, the current must in general be a 
non-monotonic function of the pumping frequency when 
there is no bias. 

It would be useful to generalize our results to the 
case of one or more strong impurity potentials, or weak 
tunnelings between two Tomonaga-Luttinger liquids; the 
technique of bosonization can be used in such situations 
also. 
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APPENDIX A: SOME MATHEMATICAL 
FORMULAE 

We need to evaluate integrals of the form 



dr 



exp(±zilT) 
((T±ia)^ -x^)^'- 



(Al) 



where fl — u!± and x — Xrp/v. Eq. (|Aip can be written 
as the sum of integrals running from to a; and from x to 
oo. We find that there are several integrals from to x 
which cancel each other. One is then left with integrals 
running from a; to oo in which one can take the limit 
a ^ in the denominator. We then use the following 
results 1591 



' ,_ Bin(»r) ^ r^^V^'^'^rn K\ 

X Jk-i/2{^x), 

COs(ilT) -y/TT f2x\^^^ ^ 



dr 



(t2 - x^)K 2 \n 



m-K) 

(A2) 



where J and Y are Bessel functions of the first and sec- 
ond kind respectively. The above equations are valid for 
x, 17 > 0, and Q < K < 1. We then use the power series 
expansion [60|] 



■J-^^> - \2) ^„ ^ > n\ r(n + :.+ l)' 



(A3) 
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and the relation 

Y,{z) = . ] , [cos{^v) Mz) - J_,(z)] (A4) 
sm yKv ) 

which is vaUd for non-integer values of v. 

Finally, the following identities involving Gamma func- 
tions are useful [601, 



r(i-z) r(^) 



sin(7rz) ' 



r(z)r(z + i/2) = ^r(2z). (A5) 
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